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Abstract: We study the evolution of the universe which contains a multiple number 
of non-relativistic scalar fields decaying into both radiation and pressureless matter. We 
present a powerful analytic formalism to calculate the matter and radiation curvature 
perturbations, and find that our analytic estimates agree with full numerical results within 
an error of less than one percent. Also we discuss the isocurvature perturbation between 
matter and radiation components, which may be detected by near future cosmological 
observations, and point out that it crucially depends on the branching ratio of the decay 
rate of the scalar fields and that it is hard to make any model independent predictions. 
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1. Introduction 

Nowadays it is widely accepted that the primordial density perturbations arc the origin of 
the temperature anisotropy in the cosmic microwave background (CMB) and large scale 
structure in the present observable universe. Various cosmological observations indicate 
that these perturbations are adiabatic and Gaussian, with an almost scale invariant spec- 
trum [1]. Interestingly, these observational facts are consistent with an earlier inflationary 
era [2]: during inflation, quantum fluctuations of a slowly rolling scalar field which domi- 
nates the energy density, the infiaton, are stretched and become classical perturbations due 
to the quasi exponential expansion of the universe. A particularly convenient quantity to 
study these perturbations is the curvature perturbation ^ on uniform density hypersurfaces, 
developed in [3], or TZc on comoving hypersurfaces which is equivalent to on large scales. 
For single field inflation cases C is known to be conserved on large scales since perturbations 
are purely adiabatic, and one can obtain the power spectrum of these perturbations with 
good enough accuracy [4]. Note that, in multi-field inflationary models, in contrast, there 
exists in general a non-adiabatic pressure perturbation and this makes ( no more conserved 
on large scales* [6]. 

*This is why the power spectrum of primordial perturbations is evaluated only after the possible trajec- 
tories of the inflaton fields coalesce in the so-called 5N formalism [5] . 



- 1 - 



In conventional inflationary models, the inflaton field is assumed to play two roles at the 
same time: it dominates the energy density during inflation and makes the universe expand 
enough to solve many cosmological problems such as homogeneity, isotropy and flatness 
of the observable universe. Also, its vacuum fluctuations are relevant for the curvature 
perturbation ( and thus responsible for the primordial density perturbations. Generally, 
the latter requirement introduce extra fine tuning into the model: for example, in the 
simplest chaotic inflation model with V{(t>) = m?(j)'^ /2, the inflaton mass m can be as large 
as O(mpi), where mpi = (SvrG)^^/^ 2.4 x 10^^ GeV is the reduced Planck mass, when 
we do not mind perturbations and try to solve other problems. However, to match the 
observed amplitude of density perturbations on large scales, we need m ~ O (lO~^) mpi, 
i.e. we need a relative fine tuning of one part over lO'' [7]. However, during inflation, any 
scalar fields with their masses being smaller than the Hubble scale acquire almost scale 
invariant fluctuations. Such fields, depending on the post-inflationary evolution of the 
universe, may later generate primordial density perturbations by transferring their almost 
scale invariant isocurvature perturbations to the curvature perturbation. 

If this is the case, i.e. in the so-called curvaton scenario''^ [9], such a field, dubbed the 
"curvaton", should satisfy several requirements: firstly, its effective mass must be light, 
i.e. less than the Hubble parameter during inflation, to produce an almost fiat spectrum 
of fluctuations and to remain sub-dominant during inflation. It should also couple very 
weakly to other fields so that its potential in the early universe is not modified appreciably. 
It is also demanded that it keeps some level of non-zero value [10] and has not yet relaxed 
to its vacuum expectation value. This is necessary to generate the appropriate amplitude of 
perturbations. These conditions are basically what the conventional inflaton field should 
satisfy, which is assumed to be responsible for the primordial density perturbations, as 
well as the enough expansion of the universe. Thus, the curvaton scenario may find its 
natural accommodation in the context of multi- field inflation [11]: for example, in a recently 
proposed scenario [12] where a number of string axion fields drive inflation, it is known 
[13] that there are a number of fields which have not yet relaxed to their minima of the 
effective potential, with their mass being very small relative to the Hubble parameter during 
inflation due to the assisted inflation mechanism [14]. 

Therefore, it is natural to consider the case where multiple curvaton fields are respon- 
sible for the generation of the curvature perturbation after inflation. The fluctuations of 
these curvaton fields are non-adiabatic in nature and thus, as mentioned above, the curva- 
ture perturbation ( does not remain constant but evolves according to the energy transfer 
between different components which constitute the universe. In this paper we study this 
general curvaton model. This paper is outlined as follows. In Sec. 2, we introduce the 
coupled equations which determine the evolution of the universe. In Sec. 3 we solve these 
equations analytically using the so-called sudden decay approximation, using a novel and 
model independent method. In Sec. 4 we apply our results of the previous sections to 
several examples and compare the analytic estimates with numerical calculations. Finally 
in Sec. 5 we summarise and present our conclusions. 

t There have been some studies on similar scenario using the decay of neutrino dark matter particles [8]. 
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2. Background equations and perturbations 



In this section we will summarise the evolution of the background quantities in a flat 
universe and show the evolution equations of the curvature perturbations of the components 
in the system of multiple curvatons decaying into radiation and matter. We assume that 
the universe is initially dominated by radiation due to the decay of the inflaton field(s) 
after inflation. 

We assume that the curvatons (crj) decay into both radiation (7) and non-relativistic 
matter (m) with constant decay rates T^^ and respectively, which are fixed by under- 
lying physics. The energy transfer equations between components are given hy^ 

Q, = -(r« + rW)p,^-r«p,, (2.1) 

Q^ = ^r»p, = ^Q^i, (2.2) 

i i 
i i 

where we have introduced the total decay width of Uj, F^*-* = + Fm , and the energy 
transfer to radiation (matter) by the decay of (Tj, Q^i {Qmi)-, respectively. Note that they 
obey the constraint of energy conservation 

^^Qi + Q-i + Qm = . (2.4) 

i 

Thus from the general continuity equation of each component including energy transfer [16], 

Pa = -^H{pa + Pa) + Qa , (2-5) 

we find that for each component 

Pi = -{3H + T(^)pi, (2.6) 
Py = -4Hp^ + Y,^'^P^, (2.7) 

i 

Pm = -SHpm + J2T^^Pi. (2.8) 

i 

Note that we can obtain the continuity equation of the total energy density by summing 
over that of each component, 

p = -3H{p + p) 

" (2.9) 



= -H ^pj + Spjn + 3 ^ Pi^ , 



*One can add the effect of dark matter freeze-out and annihilation[15], but the quahtative evolution is 
not too different. 
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where the total density p and pressure p are given by 

P = 9-1 + Pm + ^Pi, (2.10) 

i 

P = Pl+Pm + ^Pi, (2.11) 

i 

respectively. In the above we take = p-y/S and pm = Pi = 0, i.e. the equation of state of 
the curvaton fields are effectively equivalent to that of pressureless matter. 

By adopting the density parameters Cl^, fi^ and 0^, we can rewrite Eqs. (2.6), (2.7) 
and (2.8) in more convenient dimensionless forms for numerical calculation. Prom the 
Priedmann equation 



1 



the density parameters satisfy the relation 

i 

Then, Eqs. (2.6), (2.7) and (2.8) can be written as 

n'i = Qi (n^ - i?-^rW) , (2.14) 
n'^ = H-^^QiTi;^ (2.15) 

i 

n'^ = H-'J2^iT(^+^'f^rn, (2.16) 
i 

and Eq. (2.12) as 

H' = -^-^H, (2.17) 
where a prime denotes a derivative with respect to the number of e-folds, 

N = J Hdt. (2.18) 

The total curvature perturbation on uniform curvature hypersurfaces is given by 

C = -H^-^ , (2.19) 
P 

which can be written as a weighted sum of the curvature perturbation of the component 
a on the corresponding uniform density hypersurfaces [6] , 

C = E (2-20) 

a ^ 

where 

Ca = -H^-^ . (2.21) 

Pa 
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The difference between any two components gives an isocurvature perturbation [17] 



Sap = 3(Ca - C/3) • (2.22) 

The total curvature perturbation on large scales evolves as [6] 

C = ^Sp^^ , (2.23) 

P + P 

where the non-adiabatic pressure perturbation is given by 

(2.24) 

P 

Therefore, as mentioned before, ( remains constant on large scales when the perturbations 
are purely adiabatic. Prom Eqs. (2.21) and (2.23), and using the perturbed continuity 
equations of each component [17, 18], we can find that the curvature perturbations of the 
components evolve on large scales as 

_ g-^r(^)(3 + n,) _ 
2(3 + i?-ir«) ' ^ ' 



5^J7,(3 + ii-irW)ii-irW (O-C, 
-^^^^^^^(3 + J^.)(C-C.) 



(2.26) 



C= ' 3^^^ ^ (C-Ct)- (2-27) 



Here we do not solve the evolution of Cm directly, though it is straightforward to write the 
evolution equation of (m- rather, from Eq. (2.20), is calculated as 

Cm = . (2.28) 

Pm 

The reason is the existence of singularity in (rn, because there exists some moment Pm = 
when the dilution of matter due to the expansion of the universe is balanced with the 
creation of matter due to the curvaton decay^ [17, 18]. 

We may solve Eqs. (2.14)-(2.17) and (2.25)-(2.27) numerically, which would be the 
simplest way to study the cvohition of the curvature perturbation. However, we can obtain 
further insights by implementing analytic analysis. In the following section we will find the 
final curvature perturbations under the so-called sudden decay approximation [19]. 



^In fact this is the same for radiation component. However, as long as we assume that the density of 
radiation is initially high so that the universe is radiation dominated, p-y < always. 
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3. Analytic approximation 



In this section, we study the curvature perturbations under the assumption that there is 
no interaction between components until the curvaton fields decay and that the decay of 
each curvaton is instantaneous. Under this "sudden decay approximation", we can derive 
analytic estimates for the curvature perturbations associated with matter and radiation 
after all the curvaton fields decay, as we will see in this section. Note that from Eqs. (2.1), 
(2.2) and (2.3), after all the curvatons decay, there is no energy transfer between matter 
and radiation and hence (m and (j are constant, though ( will still evolve on large scales. In 
this sense, we will call these and after the decay of the curvatons as "final" curvature 
perturbations, and denote by the superscript (out). 

For our purpose in this section, we decompose the radiation and matter density ac- 
cording to the source of generation, 

P-y = P7O + y~lP7i' 

^ (3-1) 

Pm — PmO "I" / ^ Pmi ) 

i 

where p^o (Pmo) is the energy density of radiation (matter) which is due to the decay of 
the inflaton field (s) and independent of the curvaton decay, and pji (pmi) is the radiation 
(matter) density generated from the decay of ai [20]. Then, Eq. (2.10) can be written as 

P = P70 + PmO + {p^i ^fk^^^^ (^^'"' + fS)^^ 

= P7O + PmO + 'Yp'ri + 'Y Pmi , (3.2) 

i i 

where we have introduced two composite densities p^i and pmi which will play the central 
role in the discussions below. 

3.1 Matter curvature perturbation 

Prom Eqs. (2.1), (2.3) and (3.2), we can see that for the composite density Pmi, 

Qmi ~ Qmi ~^ p(i) ~ ^ ' (3-3) 

i.e. the energy transfer is zero. Moreover, since the corresponding equation of state is that 
of pressureless matter, we can write 

Pmi = -^Hpmi , (3.4) 
and therefore the associated curvature perturbation [18], 

7 TT^Pmi 5p 

mi /„ r.\ 

Qmi = -H— = ^r:^ , (3.5) 

Pmi ^P'^'' 
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is conserved on large scales. Well before the curvaton fjj decays pn 



5pr, 



so 



Cmi = Ci^"^^ ) meanwhile after cTj decays pi is negligible and thus = Cmi"*^ • Therefore we 
have 



>(out) ^ >(in) 
'mi 'i 



(3.6) 



Thus, from Eqs. (2.21), (3.1) and (3.6), we find that the final matter curvature perturbation 
after all the curvatons decay is given by 



/-(out) _ / _(in) ^(in) \ ^ (in) .(in) 

Sm ~ I A'mO SmO -p(j) Pi Si 



(in) Tm (in) 

PmO + 2^ f(j)P3 



(in) 
i ' 



(3.7) 



i=0 



where ^o"^^ = ^rnO ■ "^^^ transfer coefficient Si we have introduced above is given by 



«0 - i^mO 



(in) 
mO 



.(i) 

.7 



r-- /■ ^ 



.0) 



^(in) V- im j^(in) 



(z = 1,2,--- ,n) 



(3.8) 



So we can see that the final matter curvature perturbation is completely determined by 
the decay rate and the initial energy density p^'"^ , or equivalently, initial density parameter 
ri^'"^ of each curvaton field and that of pre-existing matter as shown above. 

3.2 Radiation curvature perturbation 

In the previous section, we could use the conservation of the curvature perturbation of 
each composite density pmi to find out the final matter curvature perturbation. This is 
possible since every p^i has no energy transfer and in addition unique equation of state. 
One may hope that similar argument is applied to the other composite density we have 
introduced in Eq. (3.2), p^j, but this is not the case. Nevertheless, p^j turns out to be an 
useful quantity to calculate the final radiation curvature perturbation as we will see shortly. 
In this section, we assume that the decay rates of the curvaton fields are different so that 
they do not decay at the same time: rather, they decay successively due to different decay 
rates. Without loss of generality we put the order of curvatons by the decay rate of each 
curvaton to satisfy T^*) > r(*+^). 

First we consider a limited time interval around the decay of the curvaton field ai, 
which is assumed to have the largest decay rate. We write a combined density of radiation 
and the curvaton field 



P7O + P-yl 



(3.9) 
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Note that although the energy transfer of is zero, its equation of state is not unique and 
thus the corresponding curvature perturbation (^^^-' evolves on large scales. Therefore, as 
mentioned above, unlike pmi we cannot simply connect the initial curvature perturbations 
in the curvaton fields to the final one in radiation, but rather we have to get through the 
moments of decay. Now we assume that until (Ji decays instantaneously there is no energy 
transfer between the curvaton and radiation. Then, p^^ before and after ai decays, which 
we write respectively 

^(1) 

P^')|before = P^O<' + ^^0<^ (3.10) 
/3^'^|after = P70^'+PTl , (3-11) 

where the superscript < 1 (> 1) means that it is evaluated before (after) ai decays, and 
these densities have the same value at the moment of decay. Since p-yi is generated only 
after cri decays, the value of p^i at the moment of decay corresponds to its initial value 
and thus 

Using the fact that both p-yo and p-yi scale as a""^, we can write the ratio p^i/p^o at late 
times, which is constant after ai decays, as 

P^:^r\-''-'/-)' r? p^^ 

The individual curvature perturbations (^q and Ci remain constant on large scales before 
o"i decays. Then, the combined curvature perturbation (^!y^^ corresponding to p^^ is written 
as 

4^)«(l-/i)C^o + /iCi, (3.14) 

where 

p(i) 

3j;^Pl /o-,r.\ 

h = p(T^ • (3.15) 

4P70 + 3j;^pi 

Here /i, the weight of C,i, solely describes the evolution of C-y^^ on large scales. After the 
curvaton cri decays, the energy density p^^ is identical to p^ at that time and has a unique 
equation of state. Hence after the decay of a\, C,^^ becomes constant on large scales until 
the curvaton with the next largest decay width begins to decay, i.e. [17] 

C7^^ = Idee « (1 - A^""'^) C,^o^ + ft'^k?^ , (3.16) 

where, using Eq. (3.13), f^^'^'^^ is given by 

»(dec) _ 3P71/P7O /„ l-7\ 

^1 -4 + 3p^i/p,o- ^ ' 
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We can take the same step for the successive curvaton decays: e.g. for o"2 which has 
the next largest decay width, we just replace 

^ I after =^ P70(new) i 

C7 ^ ^ C70(new) 1 (3.18) 

and so on. In general after i-th curvaton Oi decays, the curvature perturbation in the 
radiation component is constant until the decay of (i + l)-th curvaton, and is written as 

e,' - (1 - ff^^) er' + ff^^P , (3.19) 

where 

y(dec) ^ 'iP'ii/P-iQ ^2 20) 

4 Z)fc=0 P-yk/p^O + 3p7i/ p-fO 

Therefore, after all the n curvatons decay, we find the final curvature perturbations in 
radiation as 

^(out) ^ _ ^(dec)^ ^>n-l ^ /(dec)^(in) 

= (1 - ii^-)) (1 - /f? ) c>--^ + (1 - /(d-)) /iifcS + fk^'^'Ht^ 

n 

-E-^cf^ (3-21) 

where ^q"^^ = C,^^^ and fl^^^^ = 1. The transfer coefficient is given by 

fc=i+l \ k=i+l ) 

Tn = 4''^'=) , (3.22) 

and is completely determined once we find the ratio p^j/p^o- 
3.3 Ratio of radiation after curvaton decay 

We found in the previous section that the final radiation curvature perturbation depends 
on the ratio of the radiation generated from curvaton decay with respect to the original 
radiation component. In this section, we present a general and simple way to calculate this 
ratio analytically. 

From Eqs. (2.6), (2.7) and (2.8), we can write the continuity equations of the compo- 
nents used in Eq. (3.2) as 

Pjo = -^Hp^Q , (3.23) 
Pmo = -^Hpmo , (3.24) 

p(i) 

p^i = -AHp^, - SH^p, , (3.25) 
Pmi = -^Hpmi - SH^pi . (3.26) 
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We can solve these equation analytically and the solutions are given by 

4 



(in) 
P7O - P7O 



(in) 

PmO - PmO 



a 



(in) 



(in) / a 



(in) 



exp 



-r«(t-to) 



m (in) 



f 

J to 



exp 



dt', 



r(0 



p. 



1 — exp 



-r«(i-io)]} , 



where we have set the initial time to be to- Now, introducing [21] 

Xi = ^ t , 



z 



; i 1 1 1 1 using Eq. (3.2), the Friedmann equation, 

1 



Prai ~l" Pi 



(3.27) 
(3.28) 
(3.29) 
(3.30) 
(3.31) 

(3.32) 

(3.33) 



becomes 

i\ 2 

z 
z 

where 



X 



H 



r,(in) 







(in) ■ 
7O 



+ E 



-.(i) o(in) 



"7O 



r(^) 



70 



(in) 



-1/2 r(^) 



70 

(3.34) 



(3.35) 



TO ii-{in) ' 

r(^)t' and a prime denotes a derivative with respect to xi. We choose x\ for conve- 
nience since the dependence on this particular choice of xi is absorbed into the definition 
of xh as shown above. Finally, introducing a new variable 



_ 1/2 
y ^x^ z, 

we finally obtain the dimensionless Friedmann equation 

.'\ 2 



-3 



Jxo 



(3.36) 



(3.37) 



where the coefficients aj, Pi and /3o are defined by 



_ -1-7 -1/2 



70 



r(0 j^(in 



^0 ^ ^W,;7V2 



7O 



(in) ' 
7O 



(3.38) 
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Figure 1: Plot of the integral J^^ y exp{—Ui)dui versus x = Xi. The value changes quickly only 
near the moment of decay time, f ~ F"^, and becomes constant afterwards. 



respectively. Then the ratio of radiations which determines the transfer coefficient r, is 
given by 

^ = r ye-^^dui , (3.39) 

PlO Jxo 

where the integrand is suppressed exponentially after the curvaton decays, and the 
integral becomes almost constant. In Fig. 1, we plot this integral as a function of x = Tt. 
A large change occurs only around the decay time (x ~ 1) and soon becomes constant. We 
can see that the most significant contribution of this integral comes from the epoch around 
the moment of decay. 



3.4 Final curvature and isocurvature perturbations 

After all the curvaton fields decay, i.e. fij = 0, we are left with the overall curvature 
perturbation given by, from Eq. (2.20), 



A _ Pt A(out) I Pm ^ (out) 
40^ + SQm 



A (out) 



(3.40) 



The final matter and radiation curvature perturbations are constant on large scales and 
given by Eqs. (3.7) and (3.21), respectively. Their transfer coefficients are determined by 
Eqs. (3.8), (3.22) and (3.39). Thus, the isocurvature perturbation between matter and 
radiation components = 3(Cm — C7)) which is fixed after all the curvaton fields decay 
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so that C-y and (^m become constants, is written as 

c(out) ^of Aout) _ /■(out)\ 

= 3j2isi - nXt^ . (3.41) 

i 

A particularly simple case is when all the decay rates are the same: then, from Eq. (3.8), 
the transfer coefficient of matter curvature perturbation becomes simply 

Si = = , 3.42 

where we have assumed that initially there is no matter component. As can be seen clearly, 
the most significant contribution to the final matter curvature perturbation comes from 
the curvaton field which initially occupies the largest energy density among the curvatons. 
For ri, we only need to consider a single moment of decay since the curvaton fields decay 
at the same time. Thus, from Eqs. (3.13) and (3.15), we simply have 

(dec)^ 3^/p^^ (3.43) 
and the final radiation curvature perturbation becomes, from Eq. (3.16), 

Now, from Eqs. (3.38) and (3.39), we can see that the ratio p^i/p^o is proportional to cnj, 
which is again proportional to ^\"^\ since [integral] = J^J y ex.p[—Ui)dui will have the same 
value as discussed in the previous section. Hence, 

_ .(dec) _ 3 [integral] «^ _ 3C [integral] of"^ 

~ ' ~ 4 + 3[integral] E,- " 4 + 3C[integral] ' 

where C = T^H^'^^'^^ / ( r3/2of ^'^^^ is the common coefficient of proportionality of aj to 



nf'^K Thus, with one further assumption that the initial radiation curvature perturbation 
is negligible, i.e. C^J^'* 0, the final isocurvature perturbation is, from Eq. (3.41), 



3C [integral] 



1 - n^'^^ 4 + 3C[integral] (l - 



^f^Ct^ , (3.46) 



and thus the transfer from the initial curvature perturbation <^^^^^ is proportional to the 
corresponding initial density fraction 

4. Applications 

In this section, we apply our analytic estimates obtained in the previous section to several 
examples and compare with numerical results. 
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4.1 Single curvaton 

First we consider a simple example where a single curvaton field decays into radiation and 
matter with decay rates and , respectively. If we assume that the initial curvature 
perturbation in radiation is negligible, which is usually taken as the initial condition for 
the curvaton scenario, the radiation curvature perturbation after curvaton decay is purely 
due to the decay of the curvaton field and from Eqs. (3.21) is given by 



^(out)^^{dcc)^(i„)^ (4.1) 



»7 

where 

/.(dec) 

•'^ ~ 4 + 3p^i/p^o ■ ^ ' ' 

As discussed in the previous sections this is constant after the curvaton decay, and is 
completely determined once we find the ratio Pji/pjo- This ratio is given by Eq. (3.39) as 

^ = ai r y{ui)e-^^dui , (4.3) 

and depends only on xh and ai given by Eqs. (3.35) and (3.38), respectively. 
If initially radiation dominates, i.e. fi^™^ 1, we find that 

r(i) 

where a\ becomes identical with p of Ref. [18] in the limit V^"^ S> . In this case since 
the universe is dominated by radiation component, a oc 

and H = {2t)-^. That is, for 

small xi the solution of Eq. (3.37) is given by 

2/(xi) = (2xi)V2, (4.6) 

with = X}j/2, and we can see that y{xi) is independent of Xh [21]. Thus the curvature 
perturbation depends only on ai, which is shown by using the phase space plot in Refs. [17, 
18]. 

Furthermore, in the case that the curvaton does not dominate the density during the 
evolution, we can further approximate Eq. (4.3) analytically. From the sudden decay ap- 
proximation, we can see that ai + /3i Jl["^^(r(^)/i7("^))~^/^ ^ 1 guarantees the radiation 
domination [17], thus using y(t) ~ x^Jj^ {t/tQ)^/"^ , we obtain 



.1/2 rxi 

V2ai f ' ul^^e-'^'dui, (4.7) 
Jxo 
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where we have used Eqs. (4.4) and (4.5) and to = 1/ (2i7'^™)) for the last equality. Now let 
us take a look at the integral: it is integrated from the initial time to some later time after 
the curvaton field decays. Since we are free to choose the initial time and the integrand is 
suppressed at later times after the curvaton decay, we can take the range of integration from 
zero to infinity without loss of generality. Then the integral becomes just r(3/2) = \/7r/2- 
Hence, 



P7I 
P70 



ai ^ 1.25331q;i 



Therefore from Eqs. (4.1), (4.2) and (4.8), the final curvature perturbation 
curvaton decays is 



(out) 



A(out) 

S7 



TT (in) 



0.939986 aiQ 



(in) 



(4.8) 



after the 



(4.9) 



which is in good agreement with Ref. [18]. 

In the opposite limit (ai + /3i ^ 1), i.e. the curvaton field completely dominates the 
energy density of the universe before it decays, mostly the region of integration is the matter 
dominated epoch, thus a oc ^2/3. We can find the time of the transition from radiation 
dominated to curvaton dominated era {iljo = i^i) from sudden decay approximation, 

ttr 



2i7(i"^) 



Using this , the integral becomes 

1/2 



P7I 
P-rO 

where xtr 
we find 



1/2 
aiXj^ 



Jxo \toJ 



(in) - 
(in) 



(4.10) 



1/2 



XI 



Xtr 



2/3 



-ui 



dui 



(4.11) 



Tttr- Ignoring the contribution from the transient radiation dominated era. 



Pi'^ ^ 1/2 
P-yO 



22/3 




-Ul 



du-[ 
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(4.12) 



where we assume that initially radiation dominates the universe so that xh 
and to ~ l/(2iJ('°)). 

For the final matter curvature perturbation, assuming that initially there is no matter 
component, it is independent of the curvaton domination and is simply given from Eqs. (3.7) 
and (3.8) by 
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Figure 2: The evolution of density parameters (upper row) and curvature perturbations (lower 
row) for three cases of two curvaton decays: as shown, the energy densities of the two curvaton 
fields are sub-dominant (left panel), dominant (right panel) and comparable (middle panel) to the 
radiation energy density. The details are given in Table 1. 

i.e. it is just the same as the initial curvature perturbation in the curvaton, as shown in 
Ref. [18]. 

4.2 Two curvatons 

In this section we consider the next simplest case where there are two curvaton fields 
decaying into both radiation and matter. If we assume again that the initial curvature 
perturbation in radiation is negligible, the final curvature perturbation in radiation is, 
from Eq. (3.21), 

where 

^(dec) _ 3p^i/p^o 

^ 4 + 3/9^i//)^o ' 

J. (dec) 

^ 4(1 + p^l/p^o) +3Pt,2//070 ' 

The final curvature perturbation in matter is given by, from Eq. (3.7), 

Therefore the final isocurvature perturbation between radiation and matter is now com- 
pletely determined from Eq. (3.41). In Fig. 2, we show some examples where two curvaton 
fields decay into radiation and matter. 
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analytic 


0.883149 


0.933227 


0.889327 


numerical 


0.883149 


0.933228 


0.889328 



Table 1: The analytic and numerical results for the cases shown in Fig. 2. In the upper half of 
the table, we give the initial values used in the calculation and in the lower half we compare the 
analytic estimates with the numerical results. For analytic approximation we first solved Eq. (3.37) 
to find the density ratio Eq. (3.39) and then used Eqs. (3.20), (3.21) and (3.22). For analytic limit 
we use Eqs. (4.8) and (4.12) to find the radiation ratio in the both limits where the curvaton fields 
remain sub-dominant / dominant. For final matter curvature perturbation we have used Eqs. (3.7) 
and (3.8) for analytic estimate. Note that in the middle panel we did not use any of the analytic 
limits, since in this case the curvatons fields occupy an amount of energy density comparable to 
radiation and does not correspond to any of the limiting cases. 

If the energy density of the curvaton fields remains sub-dominant throughout the 
evolution of the universe, which would be guaranteed by the conditions 

cti + /3i < 1, 02 + /32 < 1, (4.17) 

then 

J. (dec) 

}{ « CRUl , 

ft^'^^cna^. (4.18) 

where cr = 3^^2/4 « 0.939986. 

Now it is clear that the transfer coefficients arc proportional to the initial density 
parameter of the corresponding curvaton fields. Thus with sub-dominant curvatons the 
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(4.19) 



Figure 3: The same as Fig. 2 but with five curvaton fields. The details of the parameters are given 
in Table 2. 

isocurvature perturbation is 

i 

Pi \ Jin) , ( h \ Jin) 

For the curvaton dominated case before they decay, we can take similar steps as ra- 
diation dominated one. For example in the case of the right panel of Fig. 2, the transfer 
coefficients of the matter curvature perturbation Sj are 

^1(2) = ^^, (4.20) 

where we have used V^^ /T^^'^ = F^^Vr^^^- Since the two curvatons dominate at the same 
epoch, we can use the same normalisation for y-function, thus /i and /2 are easily approx- 
imated as 

/i-i, /2-^^ = 7rxcr' (4-21) 

where in the last equality we have used /T^^^ = T^'^ /T^'^\ The transfer coefficients of 
the radiation curvature perturbation are 

ri(2) = ■ (4.22) 

The isocurvature perturbation hence almost vanishes, which is shown in the right panel of 
Fig. 2. 

4.3 Multiple curvatons 

Now we consider more general case where there exist a number of curvaton fields decaying 
into radiation and matter. It is straightforward to extract the final curvature perturbations 
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Table 2: The analytic and numerical results of in Fig. 3. As in Table 1, we show the initial 
parameters in the upper half. 



either numerically by solving Eqs. (2.14) (2.17) and (2.25)-(2.27), or analytically by using 
Eqs (3.40) with Eqs. (3.7) and (3.21). Indeed, as shown in Table 2, analytic estimates give 
good approximations to the full numerical result within an error of 0.7% (5%) with analytic 
approximation (analytic limit) . However the evolution of each perturbation could be quite 
non-trivial, as shown in Fig. 3 where we have plotted several cases with five curvaton fields. 
We can read the followings: 

• The evolution of the total curvature perturbation ( depends, not surprisingly, on 
which component dominates the energy density of the universe. During the curvaton 
fields dominates the energy density before they decay, ( is the average of (iS and 
constantly evolving during this epoch, since the curvatons are decaying into radiation 
and matter. This is clearly seen in the right panel of Fig. 3. After all the curvatons 
decay, ( follows ^7 when radiation dominates before matter begins to dominate, and 
C = Cm afterwards, as shown in Eq. (3.40). 

• Cy and Cm evolve only during the curvaton fields decay and remain constant after 
curvaton fields decay since, as mentioned before, there is no energy transfer between 
radiation and matter. Especially, since matter is assumed to be produced purely 
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due to the decay of the curvatons, is greatly affected no matter the curvaton 
fields dominate the energy density before decay or not, e.g. in the left panel of 
Fig. 3 where the curvatons never contribute significantly, their impact on is large: 
when = 0, Cm is just a weighted sum of the initial curvature perturbations of the 
curvatons and the weight Sj is basically the ratio of the corresponding curvaton energy 
density multiplied by the branching ratio to matter to the total curvaton energy 
density responsible to matter density, as shown in Eq. (3.8). For i^-y, it is noticeable 
that becomes significant only when the curvaton fields occupy significant fraction of 
total energy density before they decay, as can be compared between different columns 
of Fig. 3. This is because practically the radiation is completely generated by the 
decay of curvaton fields, making the pre-existing radiation irrelevant. 

From the discussion above, one may tempted to conclude that there will be negligible 
isocurvature perturbation between matter and radiation if the curvatons dominate 
before they decay, because they are both generated due to the decay of the curvaton 
fields. This is not true when there are a number of curvaton fields: the final isocur- 
vature perturbation is dependent on the background parameters such as curvaton 
densities and decay rates. For example, in Fig. 4, the branching ratio to matter of 
the curvaton which has the largest energy density is extremely small, i.e. 



r 



(5) 



™ ^ in-16 



10~^^ . (4.23) 



Thus, although (^^n receives contribution from the decay product of the curvaton with 
large energy density and this gives a rise of Cmi this rise is never enough to catch 
up ^7 to make Sm^^^ vanishing if the branching ratio is very small as in this case. 
This is reminiscent of multi-field inflation: in multi-field inflation, there is no unique 
prediction on the isocurvature perturbation produced during inflation. The detail 
depends on the infiaton trajectory in the field space. Likewise, generally we can 
hardly make any definite prediction on the isocurvature perturbation without the 
detail. 



5. Conclusions 

In this paper, we have studied the evolution of the universe which contains a number of non- 
interacting scalar particles (the "curvatons" (T.j) decaying into radiation (7) and pressureless 
matter (rn) after inflation. We first have written the evolution equations of the background 
densities of the components pi, p^y and pm which compose the universe and of the curvature 
perturbations of corresponding component Q C7 and Cm on fiat hypersurfaces, Eqs. (2.14)- 
(2.17) and (2.25)-(2.27). These equations can be numerically solved and give the resulting 
curvature perturbations of the components, as well as the total curvature perturbation C 
given by Eq. (2.20). 

Using the sudden decay approximation, we have obtained analytic estimates of the 
final radiation and matter curvature perturbations C-y""*"* and Cm"*^ which are in good 
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Figure 4: All the parameters are the same as the right panel of Fig. 3 except the branching ratios 
to matter: here, r|™V^^'"' is given by 10"^^, 10-"•^ 10^l^ lO^^^.s g^^^ iq-ii f^j. ^^^^i curvaton, 
respectively. In this case, we have C7°"'Vcf"^ = 0.956524 (C7°"*Vcf"^ = 0.950257 with analytic 
approx. ) and Cm"*Vcf''^ 0.602384, making 5,^°"*^ not negligible. 



agreement with full numerical results. With the composite densities p^i and pmi, given by 
Eq. (3.2), we can relate and Cm"*^ to the initial curvature perturbations associated 

with the curvatons Ci'™''- the curvature perturbation (^mi is conserved and hence the final 
matter curvature perturbation Cm"*^ has very simple relation to Cj-™"* Eq. (3.7) with the 
transfer coefficient Si given by Eq. (3.8). Meanwhile, is not constant on large scales 
since the equation of state of p^i is not unique. Nevertheless, we can find that 4°"*^ IS 
written in terms of C^^'"^ as Eq. (3.21), with the transfer coefficient rj given by Eq. (3.22). 
r, is determined once the ratio p^yi/p-yo is found, and we have found a general and model 
independent result Eq. (3.39). This might be also useful to investigate non-Gaussianity of 
the primordial curvature perturbation in the multi curvaton scenario [22]. 

We have applied our results to several different cases. The analytic estimates give 
good enough fits to the full numerical results, within an error of 0(0. 1)%. When the 
curvatons dominate the energy density of the universe before they decay, the final radiation 
curvature perturbation is significantly affected by the curvature perturbations of the 

curvatons Ci, since practically radiation is generated by the decay of the curvaton fields and 
the pre-existing radiation is irrelevant. More importantly, the isocurvature perturbation 
between matter and radiation given by Eq. (3.41) depends on the detailed decay rate of the 
curvatons: for example, in the right panel of Fig. 2, 4°"*^ and Cm"*^ are of almost the same 
amplitudes and thus isocurvature perturbation is highly suppressed. However, as shown in 
Fig. 4, when the branching ratios to matter are different for different curvatons, we may 
have significant isocurvature perturbation depending on the initial values of the background 
quantities. We can determine 5^7*'* which may be detected in the CMB observations only 
when we have detailed information on the curvaton fields. 
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